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Abstract 



Q ■ Infrared reflectance spectroscopy measurements were made on four mem- 



bers of the MFe4Pi2 family of filled skutterudites, with M=La, Th, Ce and 



^ ■ U. In progressing from M=La to U the system undergoes a metal- insulator 

"iJ ■ transition. It is shown that, although the filling atom induces such dramatic 

f~^ ■ changes in the transport properties of the system, it has only a small effect on 

On ■ lattice dynamics. We discuss this property of the compounds in the context 



of their possible thermoelectric applications. 



^. 

O . I. INTRODUCTION 

O . 

k>( ; Filled skutteruditesa are a family of compounds with the general structural formula 

c^ I MT4X12 (M=alkaline earth, rare earth, actinide; T=Fe, Ru, Os; X=P, As, Sb). Previous 

studies have shown that their physical properties are in large part determined by the M 
atoms. A variety of interesting physical phenomena, including superconductivity, antiferro- 
magnetism, fluctuating valence, and a metal-insulator transition have been observed among 
the various members of the filled skutterudites^. Recently these compounds attracted re- 
newed attention due to their potential application as thermoelectric materialsQli2l. Mahan 



and SofoO theorized that optimal thermoelectric properties might be realized for materi- 
als which possess a sharp feature in the electronic density of states, which they modeled 
as a Dirac delta function. Two different approaches have been used in attempts to realize 
this situation experimentally. One is by reducing the dimensionality of certain materials, 
as in PbTe/Pbi_2;Eu2;Te multiple quantum wellsEj. The other approach, which we discuss 
here, explores electronic f-orbitals of lanthanide atoms, that are almost completely localized 
in many lanthanide-containing compounds and which give a contribution to the density of 
states which is a Lorentzian of very narrow width. In searching for materials that satisfy the 
requirements outlined by Mahan and Sofo, we studied four members of the MFe4Pi2 family 
with M=La, Th, Ce and U as filling atoms. 

II. EXPERIMENT 

Single crystals of MFe4Pi2 compounds were grown in a molten tin flux. High purity 
(99.9% or better) pieces of the M elements (La, Th, Ce or U), Fe, P and Sn were placed in 
a quartz tube in the proportions of 1:4:20:50 and sealed under 150 mm of ultra-high purity 
argon gas. The quartz tube was placed in a furnace and the temperature was elevated to 
1050-1 150°C. The quartz tube was maintained at this temperature for about a week, and 
the temperature was than decreased at the rate of 2 °C/h. When the temperature of the 
furnace reached about 600°C, the quartz tube was dropped into water, and the product 
of the reaction, consisting of several single crystals embedded in the Sn-P solution, was 
recovered. These single crystals often grew in the habit of truncated octahedra tj^ical of 
skutterudite compounds and had dimensions somewhat less than 1x1x1 mm^ (Fig.l, top 
panel). They could be isolated after etching the Sn-P matrix with a concentrated solution 
of HCl. The crystal structure of a few of these crystals was verified utilizing Gandolfi or 
Laue back-refiection x-ray cameras. 

Electrical resistivity pdci^) measurements were made using a standard four-probe tech- 
nique. The shape of the crystals prevented an accurate determination of the sample geometry 



and hence of the absolute values of pdc(T). Instead the scale of the electrical resistivity data 
was determined from the room temperature value of the resistivity inferred from the opti- 
cal measurements. For La and Th-based samples, pdc was extracted as the best fit to the 



Hagen-Rubens formula R{ijj) = 1 — J2ujpdc/T^, whereas for Ce- and U-based samples it was 
obtained by extrapolating the optical conductivity down to zero frequency. 

The electronic structure and charge dynamics in MFe4Pi2 were studied using infrared 
and optical reflectance spectroscopy. Near normal incidence reflectance R{uj) of single- 
crystal samples of MFe4Pi2 was measured at UCSD in the frequency range 50-30,000 cm~^ 
(approximately 5 meV - 3 eV). Single crystal sample of LaFe4Pi2 is shown in the top panel 
of Fig.l. In our experiments we have collected light reflected by one of the facets of the 
crystal. To obtain the absolute value of R(c(j) samples were coated in-situ with gold or 
aluminum in the optical cryostat and the spectrum of a metal-coated sample was used as 
a reference. This procedure yields reliable absolute values of R{uj) and does not require 
ambiguous corrections for diffuse reflectancet^. Filled skutterudites provide an excellent 
illustration of the importance of performing optical measurements on natural surfaces. The 
two La-based samples that were analyzed differed only with respect to the condition of their 
surface: one of the surfaces was polished, the other was not. A scanning electron microscope 
(SEM) picture of the sample with natural, unpolished surfaces is shown in Fig.l (top panel) 
together with the far infrared (FIR) reflectance of both polished and unpolished samples 
(bottom panel). Obviously, the polishing procedure significantly damages the surface layer 
of the sample, as evidenced by the following features in the reflectivity data for the polished 
sample: (1) a lower overall free-electron reflectance and (2) weaker peaks associated with 
phonons due to the lack of crystalline long-range order in the surface layer. Fortunately, 
the procedure for measuring reflectance described above allows us to work with samples of 
submillimeter size, thus eliminating the need to produce larger surfaces by polishing. The 
measured frequency interval was extended up to lOeV at ETH on polished specimens of 
LaFe4Pi2 and UFe4Pi2. 

The complex conductivity a{uj) = (Ti{uj) +i(T2{u!) and complex dielectric function e(co') = 



ei{uj) + ie2{uj) were obtained from R{u!) using a Kramers-Kronig analysis. The uncertainty 
of the cr(co') and e{uj) spectra due to both low- and high-frequency extrapolations required 
for the Kramers-Kronig analysis are negligible in the frequency range where the actual data 
exist. 

III. RESULTS 

Figure 2 shows the temperature dependence of the electrical resistivity pdc{T). The 
LaFe4Pi2 and ThFe4Pi2 compounds behave like metals, having positive temperature coef- 
ficients of resistivity. The residual resistivity ratios RRR = p{300K) / p{A.2K) for the two 
samples are 94 and 1.8, respectively. On the other hand, the resistivities of CeFe4Pi2 and 
UFe4Pi2 have negative temperature coefficients typical of insulators and semiconductors. 
The resistivities change by as much as seven orders of magnitude between room tempera- 
ture and 4.2 K. 

Presented in Fig. 3 are the raw reflectivity data at room temperature for each of the com- 
pounds. As anticipated from transport measurements, UFe4Pi2 shows insulating behavior. 
The absence of the free electron background allows detailed analysis of the lattice dynamics. 
Another noticeable feature is relatively high reflectivity (~ 40 %) in the mid-infrared region 
(500 - 5,000 cm~^). In most insulating compounds R(co') does not exceed 10 % in this range. 
We propose that this is a result of the high frequency contribution to the dielectric constant, 
eoo- This behavior seems to be common for all insulating skutterudites, and will be discussed 
further below. 

For CeFe4Pi2, the other known nonmetallic member of the skutterudite family, the re- 
flectances of the two samples analyzed were the same at frequencies u >300 cm~^, but were 
quite different from each other in the far infrared region. Sample No.l behaved like a semi- 
conductor (some free electron background is present), whereas sample No. 2 was insulating. 
This sample-to-sample variation was already evident in measurements of both transport and 
thermodynamic quantities^. Since our optical experiments probe the bulk of the crystalli^ 



and results show significant variations between samples, we conclude that this is a bulk 
effect. Similar to UFe4Pi2, this material exhibits a high reflectance (~ 50 %) in the mid- 
infrared part of the spectrum. The increase in reflectivity in the near-infrared and visible is 
suggestive of interband transitions. 

The third member of the family, ThFe4Pi2, displays metallic behavior. The free-electron 
reflectance is much higher here than in CeFe4Pi2 and UFe4Pi2, but the relatively low plasma 
edge (around 3,500cm~^) is an indication of a low carrier concentration. The increase of 
the reflectivity in the near-infrared and visible parts of the spectrum is again indicative of 
interband transitions. 

Finally, LaFe4Pi2 is also metallic and exhibits a far- infrared reflectivity that decreases as 
temperature increases (data not presented in this paper). From transport measurements it 
is known that this compound is a superconductor with a transition temperature of around 
4.2 KB. The general shape of the reflectivity is the same as in ThFe4Pi2. The only difference 
between the two is in the position of the "plasma minimum" in R(ci;), which is shifted to 
about 4,000 cm~^ in LaFe4Pi2 (see Fig. 3). 

The next step in the analysis is to perform a Kramers-Kronig transformation on the raw 
reflectivity data. To do that one needs to extrapolate the reflectivity to the regions where the 
actual data are not available. For electrically conducting samples (LaFe4Pi2 and ThFe4Pi2), 
high FIR reflectance suggests using a Hagen- Rubens (HR) extrapolation (commonly used for 



metals): R{uj) = 1 — J2iupiic/'Tc, where pdc is the only parameter of the model. The natural 
way to proceed is to use pdc from transport measurements and calculate R{uj). But since pdc 
was not available, we inverted the procedure: we found the best fit to the reflectance data and 
used it to calculate the absolute values of our transport measurements (Fig. 2). The best fits 
for LaFe4Pi2 and ThFe4Pi2 were 85 and 200 pQcm, respectively. A power law extrapolation 
R{uj)r^l/uj^ was used for high frequencies. For the CeFe4Pi2 and UFe4Pi2 samples, we 
extrapolated the low frequency part with a straight line: for the semiconducting CeFe4Pi2 
sample perfect reflectance at zero frequency was assumed (suggested by the tendency of the 
data), whereas for the insulating CeFe4Pi2 and UFe4Pi2 compounds we assumed a constant 



reflectance (horizontal straight hne). For the high frequency part the same power law as 
before was used. 

Using these extrapolations, the Kramers-Kronig transformation was performed and the 
frequency dependence for the complex conductivity cr{u!) and complex dielectric function 
e{uj) were obtained. Figure 4 shows results for the real part of the conductivity ai{uj). 
Characteristics anticipated from reflectance data are now obvious: LaFe4Pi2 and ThFe4Pi2 
are metallic, one of the CeFe4Pi2 samples is semiconducting while the other is insulatingli3, 
and, flnally, UFe4Pi2 is insulating. 

Conductivity spectra for LaFe4Pi2 and ThFe4Pi2 are clearly Drude-like, but we failed 
to obtain good flts to the formula <Ji{uj) = o"o/(l + uj^t^). Instead, spectra were analyzed 
within the "extended" Drude model. Here one allows for the frequency dependence of the 
scattering rate 1/r. It can be shownllB that causality requires the effective mass to be 
frequency dependent as well. The spectra of the quasiparticle effective mass m*{u)/me 
{■rrie - electron band mass) and scattering rate 1/t{uj) can be obtained from the optical 
conductivity: 
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where the plasma frequency u^ = Anne's/me is estimated from the integration of o"i(u;) 
up to the frequency corresponding to the onset of the interband absorption. When these 
equations are used to describe the response of a conventional metal one flnds that both the 
scattering rate and effective mass are independent of frequency. The results of the analysis 
for LaFe4Pi2 and ThFe4Pi2 are presented in Fig.5 (top and bottom panel, respectively). For 
both compounds the effective mass (right panels) is only slightly enhanced {rrf/rrie <5) and 
almost temperature independent. The scattering rate (left panels) is linear in frequency (up 
to 3000 cm^^ for LaFe4Pi2 and 2000cm^^ for ThFe4Pi2) and is temperature dependent. The 
values obtained by extrapolating the curves down to zero frequency {dc limit) are plotted in 



Fig. 2 as full squares (for LaFe4Pi2) and full circles (for ThFe4Pi2). Such a good scaling is 
not surprising if one considers the relation p^c = vn* jne^r and the fact that effective mass 
m* and carrier concentration n are temperature independent. A huge factor of 50 difference 
in the residual resistivity ratio RRR between La and Th-based samples is consistent with the 
temperature dependence of the scattering rate. At low temperature the phonon contribution 
is strongly suppressed and the large residual scattering in this region for ThFe4Pi2 indicates 
that the La-based sample was much "cleaner". 

Figure 6 displays the real part of the complex dielectric function e\{uj) up to 3,000 cm~^ 
(0.3eV) for all the samples, where contributions from the interband transitions at higher 
frequencies are not shown. Also shown are the resonance plasma frequencies lo* (determined 
as zero crossings of ei(c<j)) for the LaFe4Pi2, ThFe4Pi2 and semiconducting CeFe4Pi2, and 
high-frequency contributions e^o for the insulating CeFe4Pi2 and UFe4Pi2. 

IV. DISCUSSION 

A. Electronic Properties 

Electronic properties of both filled and unfilled skutterudites have been discussed in a 
number of publication3iHi3il3. The filling atom is always assumed to be positively charged 
and only weakly bound to the transition-metal pnictogen polyanion. The bonds within the 
polyanion are considered to be covalent. A simple charge counting method is sometimes 
used to predict the properties of the compounds. According to this method, the polyanion 
[Fe4Pi2]'^~ should be isoelectronic with C0P3 (C04P12) which is semiconducting. The elec- 
tronic properties of MFe4Pi2 will therefore depend on the valence state of the filling atom 
M. Metallic behavior is indeed observed in MFe4Pi2 where M = La, Pr, Nd, Sm, or Eu 
is trivalent, and semiconducting behavior in CeFe4Pi2 and UFe4Pi2, where U and Ce may 
be tetravalent. However, there is growing consensus that f-electron hybridization effects 
better explain the semiconducting behavior in latter materialsEHiHHl, as we discuss later in 



this paper. Furthermore, the charge counting arguments fail to predict the correct behavior 
for ThFe4Pi2, in which Th is beheved to be tetravalentu and is therefore expected to make 
the compound semiconducting (electronically saturated). Both the frequency (Fig.4) and 
temperature (Fig. 2) dependence of the conductivity clearly demonstrate metallic behavior. 
This simple charge counting method thus seems insufficient to explain all the experimental 
results, and more detailed band structure calculations are necessary. 

The absolute values of the resistivity data for the metallic samples (Fig.2, top panel) are 
approximately two orders of magnitude higher than in typical metals. Apparently the curve 
for the Th-based sample exceeds the value of 150 fiQcm, which is the so-called lofFe-Regel 
limit. In conventional metals it corresponds to the regime when the carrier mean free path 
is approaching the lattice constant. The value of 150 fiQcm is calculated for conventional 
metals, where n ~ lO^^cm"^, r ~ 10~^^ s and m* /me = 1. From the plasma frequency Up we 
estimated the carrier concentration in the Th-based sample to be approximately two orders 
of magnitude lower {n ~ 10^°cm~^) which places the loffe-Regel limit for this compound 
much higher. The spectra of the dielectric constant (Fig. 6) of LaFe4Pi2 and ThFe4Pi2 also 
exhibit the characteristics of "poor" metals. The zero crossing of the dielectric constant, 
which determines the plasma resonance frequency: 



^P = -7= = \ (3) 



is in the range 0.1-0.3eV, confirming a low carrier concentration for both compounds. 

Besides having low carrier concentration, these compounds also seem to have an unusual 
frequency dependence of the carrier scattering rate: 1/r ~ u (Fig.5, left panels). We believe 
that this is caused by the interband transitions that are at relatively low frequencies (around 
10,000cm^^, Fig.4) and contribute considerably even below lOOOcm^^. The interband tran- 
sitions also affect the effective mass: they drive it below the free electron mass, and even to 
negative values (Fig.5, right panels). 

On the other hand, UFe4Pi2 and CeFe4Pi2 display all the characteristics of insulators 
or semiconductors. High values of eoo (17 for UFe4Pi2, 31 for CeFe4Pi2 No.l and 32 for 
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CeFe4Pi2 No. 2, see Fig.6) are typical for hybridization gap semiconductoraH^o, as explained 
later in the text. This results in the anomalously large mid-infrared reflectance that was 
mentioned above. Model calculations of reflectivity for several different values of eoo are 
presented in Fig. 7. The other parameters of the model, obtained as the best fit to the 
reflectance spectrum of the U-based sample, were kept constant. It is apparent that the 
increase of eoo leads to the enhancement of the absolute value of R(u;) in the mid-infrared. 
All conductivity spectra show a broad maximum centered around 1.5eV (Fig.4), a signa- 
ture of an interband transition. At this moment, band structure calculations are available 
for Ce and La-based compounds onlynHiH3, and they do not show any characteristic feature 
at this energy range. For CeFe4Pi2 a gap at the Fermi level of about 0.34eV is predicted^, 
but our spectrum indicates a somewhat smaller value of about O.lSeV. This value is also 
consistent with the activation energy obtained from fitting resistivity measurements to an 
activation conduction form, p=poexp(AE/kBT), in the temperature range 80K < T < 150K 
which yields an activation energy AE of around O.lSeVH. One possible origin of the insulat- 
ing behavior which has been previously put forwarda is hybridization between the f electron 
orbitals present in UFe4Pi2 and CeFe4Pi2 with the conduction electrons which opens a gap 
in the electronic density of states at the Fermi energy. This would be analogous to the 
hybridization gap seen in the intermediate valent systems SniBg and SmSc3, except that the 
latter materials have considerably smaller activation energies than CeFe4Pi2. Band structure 
calculationsE3 also suggest that Ce 4/ states hybridize strongly with Fe 3d and phosphorus p 
states in the vicinity of the Fermi level. The notion of hybridization is also supported by the 
experimental findings presented in this paper. Namely, hybridization gap semiconductors are 
predicted^ to have: (1) high values of eoo and (2) anomalously large values of Born effective 
charges (explained below). We indeed find experimental evidence for both of these effects. 
As La and Th do not have any f electrons, no hybridization takes place, and thus LaFe4Pi2 
and ThFe4Pi2 display metallic behavior. Clearly, additional theoretical and computational 
work is required to explain these experimental findings in metallic skutterudites. 



B. Phonon Spectra 

Filled skutterudites crystallize in the cubic structure (space group Im3 - T|c3). Applying 
the correlation methods, the following decomposition into the irreducible representations 
is obtained: F = 2Ag + 2Eg + AFg + 2Au + 2£'„ + SF^. Eight of these modes are infrared 
active (F„ species), eight are Raman active (A^, Eg and F^ species) and the other are neither 
infrared nor Raman active. All infrared active modes are triply degenerate. 

Frequencies and oscillator strengths of all phonon modes identified in the conductivity 
spectra are given in Table 1. The numerical values are obtained as the best fits in the 
classical oscillator modelcS: 

o-ph[uJ) = -j— 2^ T^ jT^-^ — 2~^ (.4) 

47r ^ [ujj — oj^Y + '^jOJ^ 

Although the group theory analysis predicts eight infrared active modes, we were able to 
identify only six from the spectra of the insulating samples. The oscillator strength of the 
other two was probably below our detection limits. Due to a high free-electron contribution 
to the conductivity, only four phonons are identified in the Th-based sample and three in 
the La-based sample. Interestingly, except for the lowest lying mode, all peaks are at nearly 
the same frequency in all compounds. We propose below that they are due to the Fe4Pi2 
sublattice. The filling ion M is believed not to be very strongly bound to this sublattice and 
is free to rattle within the "cage". Weak binding could be the reason for the high oscillator 
strength of this mode (Table. 1) The most interesting feature of the phonons is that even 
though the conductivity (carrier concentration) changes by almost three orders of magnitude 
going from the U to La samples, the oscillator strength of the phonon modes is only weakly 
affected (Table 1). Despite the fact that the conductivity of "metallic" samples is enhanced, 
the carrier density is still below 0.2 electrons/unit cell which is much lower than the effective 
oscillating charge. 
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C. Effective Charges 

The discussion in Section IV. A suggests that the metal-insulator transition in the 
MFe4Pi2 series is in accord with the opening of a hybridization gap in the electronic density 
of states. Another predicted feature of the hybridization gap semiconductors are anoma- 
lously high values of the effective dynamic chargesn3. Our calculations of the Born dynamical 
charges e^(M) for U and Ce-based samples are indeed in agreement with this prediction. 
The Born effective charge can be determined from the formula: 

T.[-LO^--TO,) = y-,^^^ (5) 

where ujlo,j and ujto,j are frequencies of the longitudinal and transverse optical phonons, 
respectively, V^ is the volume of the primitive cell, Z is the number of formula units in the 
primitive cell (Z=2) and m is the reduced mass of the primitive cell. 

In accord with our assignment of phonons in the previous section, we calculate e^ for 
the filling ion and Fe4Pi2 sublattice separately. For the former, the sum on the left-hand 
side of Eq.5 has only one term (the lowest lying phonon), whereas for the latter it runs over 
all the remaining peaks. This procedure gives the following results for the Fe4Pi2 sublattice: 
e^=8.9 in UFe4Pi2 and 11 in CeFe4Pi2. Such high values for a predominantly covalent lattice 
(covalent bonding is expected within the Fe4Pi2 polyanion) indicate that hybridization plays 
an important role. The dynamical charge associated with the filling ion is also anomalously 
high. For U and Ce-based samples we find e^ of 8.9 and 8.7, respectively. Obviously, such 
high values of the transverse dynamical charge (i.e., the one that couples to the external 
radiation and that is measured in optical experiments) cannot be due to the valence state 
of the ionically bonded filling atoms. 

D. Thermoelectric AppUcations 

There has been an enormous amount of research in the field of thermoelectric materials 
in the recent years, with emphasis on novel materials concepts. The figure of merit of a 
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material for thermoelectric applications is measured by the dimensionless quantity: 

ZT = — — T (6) 

where T is temperature, a is the electrical conductivity, S is the Seeback coefficient and Ke 
and K/ are electronic and lattice contributions to the thermal conductivity. For the novel TE 
materials to be competitive with conventional cooling systems, ZT ?^ 3 is required^. The 
ffiled skutterudites (e.g., LaFe3CoSbi2 and CeFe3CoSbi2) attracted attention as potential 
thermoelectric materials because it was shown that they exhibit large Seebeck coefficients 
typical of semiconductors yet have moderately high electrical conductivity and also very 
low lattice thermal conductivity /t^Btjc^. The latter property has been attributed to the 
localized, incoherent vibrations (rattling) of the filling ions within the atomic " cages" formed 
by the other atoms. The hypothesis was based on ultrasound, neutron diffraction, and 
specific heat measurements, and on the fact that /t/ for the unfilled skutterudite CoSba was 
several times larger at all temperaturesE3. 

We believe that the lowest lying phonon (which varies in frequency from compound to 
compound) is this rattling mode. All other phonons can be attributed to the dynamics of the 
Fe4Pi2 sublattice, which explains their similarity in all four compounds (both in frequency 
and in strength). It has been recently suggestedE^I'El that the rattling of the filling ion can 
be treated as a localized vibrational mode. We note that in UFe4Pi2 and CeFe4Pi2 the 



frequencies of the lowest lying mode scale as Jk/M, where M is the mass of the filling 
ion and k is the force constant, found to be the same in Ce and U-based crystals. This 
scaling suggest that the filling ion is weakly bound to the lattice. If the ffiling ion were 
strongly bound to the lattice, one would expect the force constant to be quite different for 
different ions. Based on the calculated value of the force constant we estimate the phonon 
frequencies for the same mode in ThFe4Pi2 and LaFe4Pi2. They are expected to be around 
82 and 141cm~^, respectively (indicated by arrows in Fig.3). The free electron background 
in raw reflectance data in this frequency range for these two compounds is too high to allow 
detection of any phonons. 
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In order to produce viable thermoelectric materials, careful tuning of the com- 
pound's carrier concentration is required. Several alloying schemes are possible within 
the filled skutterudites. Compounds of the form M(Ta;T'4_2,)Xi2, MT4(XyX'i2-j,) or even 
M(Ta;T'4_a;)(Xj,X'i2_y) havc been successfully synthesizedEMES. Recently compounds of the 
form Lai^a;Cea;Ru4Pi2 (0 < x < 1) were produced at high pressureg23, demonstrating the 
ability to vary x continuously over the entire alloying range. We believe that this may 
allow us to design a material whose electrical and thermal conductivities are optimal for 
thermoelectric applications. 

Finally, we want to discuss our experimental results from the point of view of Mahan 
and Sofo's theorytll. They point out that all transport coefficients relevant for thermoelectric 
materials {a, S and Kg) are functions of the electronic density of states Ne/(E). They find, 
in the case of a material with the minimum lattice thermal conductivity k; ^ 2mW/{anK), 
that a sharp feature (theoretically a Dirac delta function) in the electronic density of states 
located 2.4kBT away from the Fermi level is needed to maximize ZT. This is equivalent to 
having a resonance approximately 60meV (at room temperature) above or below the Fermi 
level. None of the materials we have analyzed in this paper shows any characteristic feature 
on this energy scale (Fig.4). We are currently studying another member of the skutterudites 
family, YbFe4Sbi2Q, and find interesting charge dynamics in the anticipated energy rangeE^. 

V. CONCLUSIONS 

We have presented evidence in this paper that transport and optical properties of the 
filled skutterudites are primarily determined by the filling ion. The experimental results 
support the idea that hybridization plays an essential role in determining the properties 
of the MFe4Pi2 compounds. From the optical spectrum of samples with M = Ce and U, 
we inferred the following properties: (1) a gap in the electronic density of states, (2) a 
large value of the dynamical effective charges, and (3) an enhanced dielectric constant at 
high frequencies Coo- Localized, incoherent vibrations ("rattling") of the filling ions were 
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successfully modeled as Einstein oscillators, in agreement with recent reporta^EllEj. Atomic 
substitutions at this site will not only tailor the electrical conductivity of the material but 
may also increase the phonon scattering due to the rattling of the various filling ions. We 
believe that this method (atomic substitution at the filling atom site), combined with alloying 
at the Fe and P sites, opens possibilities for significant improvements of the thermoelectric 
coefficient. Further research, especially measurements of thermal conductivity, is needed to 
verify these ideas. 
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FIGURES 
FIG. 1. An SEM picture of as-grown LaFe4Pi2 single crystal (top panel). Reflectance spectra 

(bottom panel) show the importance of natural surfaces; polishing damages the surface layer of 

the sample, creating defects that increase resistivity (lower the reflectance). The phonons are also 

much weaker due to the lack of long range order. 

FIG. 2. Temperature dependence of the dc resistivities of MFe4Pi2. The absolute values are 
obtained by extrapolating the optical conductivity down to zero frequency (for U- and Ce-based 
samples) and from the Hagen-Rubens fits to the reflectance (for Th- and La-based samples) Also 
shown is the loffe-Regel limit calculated for n~10^^ cm~^ (typical for conventional metals). In 
ThFe4Pi2 the carrier concentration is two orders of magnitude lower, placing the loffe-Regel limit 
much higher. The carrier scattering rate (right-hand axis in the top panel), is shown as full 
circles for LaFe4Pi2 and full squares for ThFe4Pi2. This values are obtained by extrapolating the 
scattering rate to zero frequency (see Fig. 5). 

FIG. 3. Reflectance spectra of four members of the skutterudite family. The arrows on the 
graphs for La and Th samples indicate the positions of the "rattling" mode, calculated from the 
simple harmonic model (see text for details). 

FIG. 4. Real part of the optical conductivity as obtained from a Kramers-Kronig analysis. All 
spectra show a broad maximum around 1.5eV, which we attribute to an interband transition. The 
low frequency part of the CeFe4Pi2 spectra (thin lines) is obtained by fitting the reflectance data. 

FIG. 5. Frequency dependent effective mass (right panels) and scattering rate (left panels) 
obtained within the "extended" Drude formalism of Eq.l and Eq.2 (see text), applied to LaFe4Pi2 
(top panels) and ThFe4Pi2 (bottom panels). Note a slight mass enhancement and linear frequency 
dependence of the scattering rate. 



FIG. 6. Real part of the dielectric function, shown over a narrower frequency range (up to 
about 0.3eV). Again, thin lines are used for the simulated region (for CeFe4Pi2)- The low values 
of the plasma resonance frequency to* are an indication of low carrier concentrations. The high 
frequency contributions to the dielectric function (coo) for insulating Ce and U-based samples are 
also shown. High values are typical of hybridization gap semiconductors. 

FIG. 7. Model calculations of the reflectance for UFe4Pi2, shown for several values of the free 
parameter eoo- The phonon parameters I'j and pj for these curves are given in Table 1. Apparently, 
an increase in Eoo leads to an increase in the mid-infrared reflectance. 
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TABLES 
TABLE L Frequency loj in cm^"'^ and reduced oscillator strength pj = Sj/Attujj (in parentheses) 
of infrared active phonon modes at q~0 obtained by a fit to the classical oscillator model. 
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